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Abstract 

Within the Thermal Wave Model framework a comparison among Wigner func- 
tion, Husimi function, and the phase-space distribution given by a particle tracking 
code is made for a particle beam travelling through a linear lens with small aberra- 
tions. The results show that the quantum-like approach seems to be very promising. 
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1 Introduction 



A satisfactory description for the transverse dynamics of a charged particle beam is a 
fundamental requirement for the new generation of accelerating machines working at very 
high luminosity. To this aim, the recently proposed Thermal Wave Model (TWM) [|l| 
could represent an interesting framework for better describing the dynamics of charge 
particle beams. 

According to this model, the beam transport is described in terms of a complex func- 
tion, the beam wave function (BWF) , whose squared modulus gives the transverse density 
profile. This function satisfies a Schrodinger-like equation in which Planck's constant is 
replaced by the transverse emittance. In this equation, the potential term accounts for 
the total interaction between the beam and the surroundings. In particular, in a particle 
accelerator the potential has to take into account both the mult ipole- like contributions, 
depending only on the machine parameters, and the collective terms, which on the con- 
trary depend on the particle distribution (self- interaction). 

In transverse dynamics, TWM has been successfully applied to a number of linear and 
non-linear problems. In particular, it seems to be capable of reproducing the main results 
of Gaussian particle beam optics (dynamics for a quadrupole-like device as well as 
of estimating the luminosity in final focusing stages of linear colliders in the presence of 
small sextupole- and octupole-like deviations 0. In addition, for the case of transverse 
dynamics in quadrupole-like devices with small sextupole and octupole deviations, TWM 
predictions have been recently compared with tracking code simulations showing a very 
satisfactory agreement 0. 

More recently, in the framework of Nelson's stochastic mechanics the Schrodinger-like 
equation of TWM has been derived 0]. 

This paper concerns with the transverse phase-space behaviour of a charged particle 
beam passing through a thin quadrupole with small sextupole and octupole deviations, in 
order to find the most suitable function to describe the phase-space distribution associated 
with the beam in the framework of TWM. 

By using the BWF found in Ref. [^, we compute the corresponding Wigner function 
(WF) 1^ which, according to the quantum mechanics formalism, should represent the 
most natural definition for the beam distribution function. Unfortunately, due to the 
uncertainty principle, the WF is not positive definite, hence for the above BWF we com- 
pute also the corresponding Husimi function (HF) 0], the so-called Q-function, which is 
alternative to WF and is in fact positive definite. A comparison of both these predictions 
with the results of a particle tracking simulation is then performed in order to select the 
most appropriate definition of the beam phase-space function in TWM. 

The paper is organized as follows. A brief presentation of BWF determination is given 
in section 2, whilst in section 3 we give the definition of both WF and HF. In section 
4, the comparison of the theoretical predictions with the simulations results is presented. 
Finally, section 5 contains conclusions and remarks. 
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2 Beam Wave Function in presence of small sextupole 
and octupole aberrations 

Let us consider a charged particle beam travelling along the 2;-axis with velocity (3c {(3 ~ 
1), and having transverse emittance e. 

We suppose that, at 2; = 0, the beam enters a focusing quadrupole-like lens of length / 
with small sextupole and octupole deviations, then it propagates in vacuo. In this region, 
if we denote with x the transverse coordinate (1-D case), the beam particles feel the 
following potential 

where ki is the quadrupole strength, k2 is the sextupole strength and ^3 is the octupole 
strength, respectively. Note that ?7(x, z) is a dimensionless energy potential, obtained 
dividing the potential energy associated with the transverse particle motion by the factor 
mo7/3^c^, where mo and 7 are the particle rest mass and the relativistic factor [1 — 
respectively. 

As already stated, in the TWM, the transverse beam dynamics is ruled by the Schrodinger- 
like equation 

• e2 92 ^ , ^ 

The z-constancy of the integral |\E'(x, 2;)^ dz, which is a consequence of the reality 
of U{x,z) in (0), suggests to interpret z)p as the transverse density profile of the 
beam. Hence, if is the total number of the beam particles, X{x,z) = N z)p is 
the transverse number density. 

We fix the initial profile as a Gaussian density distribution of r.m.s. ctq, which corre- 
sponds to the initial BWF 

Provided that aok2/{3ki) ^ 1 and aQks/{12ki) <^ 1, the quantum mechanics formalism 
for tzme-dependent perturbation theory applied to (0) for the case of a thin lens {y/k^l ^ 
1) allows us to give an approximate first-order normalized BWF in the configuration 
space. At the exit of the lens {z = I) it reads ^ 

^(xj) = ^oix) exp f -^^^ ) (1 - tSiu) Ho ( - hJ ^ ^ 




V2 



(4) 



where r = (jQA"2/6e, and oj = ctq A'3/24e, with Ki = kil {i = 1,2,3), the integrated 
aberration strengths. Remarkably, Eq. shows that, due to the aberrations, the BWF 
after passing the kick stage is a superposition of only five modes, according to the simple 
selection rules due to (|l]). Thus, after a drift of length L {L^ I) in the free space we get 
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(see also 0) 



exp 



■40-2 (L) 



exp 



;i - iSuj) Ho 



X 



X 



V2 \y/2a{L), 



X 



2v^ V V2(r{L) J ^ \ V2(x{L) ^ 



X 



i44>{L) 



where 



a{L) 
1 

R(L) 

m 



+ Kfa'^ (L - ly - 2K^Gi{L - /) + 



,i8</.(L) 



1/2 



(5) 



1 da 

o dz 



z=L 



-— { arctan 
2 I 



Ml 



+ 



+ arctan 



(6) 



Correspondingly, the Fourier transform of \E'(x, 2;) 

2;) = — / \E'(x,z) exp(— ipx/e) dx 

\ I'K J -00 



(7) 



is the BWF in the momentum space. 

Consequently, \^{x, /)p and L)p give the transverse particle density profiles ak, z = 1 
(after lens) and at z = L (after drift), respectively, whilst |$(p, /)P and |$(p, -L)p give the 
momentum distributions of the particles, at 2; = / and at ^ = L, respectively. 

In Ref. the configuration-space and momentum-space distributions are reported for 
some significative values of parameters cxo, e, i^i, K21 K^, and L. They are compared 
with the corresponding results obtained from particle tracking simulations showing a very 
satisfactory agreement. 

In the next section we extend this analysis, made separately both in configuration 
space and in momentum space, by introducing the appropriate distribution function in 
the context of TWM approach. 

Since the approximate solution @) and (H) for BWF is given for small sextupole 
and octupole deviations from a quadrupole potential (harmonic oscillator), the present 
analysis falls in the semiclassical description of the wave packet evolution (WKB theory) 
extensively treated in Ref. 0. 



3 Phase-space distributions 

According to Quantum Mechanics (QM), for a given BWF "^{xjz) we can introduce the 
density matrix p as 

p{x,y,z) = -^{x,z)^*{y,z) , (8) 
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which, in the {hra\ and \ket) Dirac's notation, is associated with the the following density 
operator 

p = \^)m ■ (9) 

Note that p has the following two properties 

i) probability conservation 

Tr(p) = l ; (10) 

ii) hermiticity 

pt = p . (11) 

On the basis of this density matrix definition, we can define the relevant phase-space 
distributions associated with the transverse beam motion within the framework of TWM. 



3.1 Wigner function 

One of the widely used phase-space representations given in QM is the one introduced 
by Weyl and Wigner. In this representation, by simply replacing Planck's constant with 
e, the phase-space beam dynamics can be described in terms of the following function, 
called Wigner function (WF) 




namely, by virtue of (H) 




It is easy to prove from (p!3| ) that 

/+00 r+oo 
/ W{x,p,z) dx dp=l , (14) 
-oo J —oo 

/ + 00 
W{x,p,z)dp , (15) 

and 

/ + 00 
W{x,p, z) dx , (16) 
-oo 

where //(p, z) represents the transverse momentum space number density (note that ri{p, z) = 
N |$(p, 2:)p). Eqs. (|T^) and (|TBp show that times W{x,p,z) is the phase-space distri- 
bution function associated with the transverse beam motion. Consequently, A and r] are 
its configuration-space and momentum-space projections, respectively. 

The averaged- quantity description can be done by introducing the following second- 
order momenta of W 

/4-00 i"-f oo 
/ W{x,p,z) dx dp= (x^) , (17) 
-oo J —oo 

/+00 r+oo ]_ 
/ xp W{x^p, z) dx dp = -{xp + px) , (18) 
-oo J —oo Z 
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/+00 r+oo 
/ p"^ W{x,p,z) dx dp= (p^) . (19) 
-oo J —oo 

They are connected with the geometrical properties of the phase-space ensemble associated 
with the beam; here we have considered the case (p) = (x) = 0. 

According to QM and simply replacing h with e, we can easily say that W{x,p, z), 
defined by (|1^) when the BWF \Ef(x, 2;) is solution of (Q) (with U arbitrary Hermitian 
potential), satisfies the following Von Neumann equation 



^ + ^ + ^ (u 

dz ^ dx e I 




W = Q . (20) 



Thus, in order to perform a phase-space analysis, one has two possibilities: either 
to solve @j for ^>{x,z) and then, by means of the Wigner transform (0), to obtain 
W{x,p,z), or to solve directly (pO|) for W{x,p, z). 

Although for the potential given by (P the exact solution of the Schrodinger-like 
equation (^ is unknown, time-dependent perturbation theory can be applied to give, at 
any order of the perturbative expansion, the approximate solution : for example, (^ is 
the first-order approximate solution of @ for the potential ([l|). 



Also for the Von Neumann equation, (pO]) , the exact solution is not available in the case of 
potential (|l|), with the exception of the case with k2 = k^, = (pure quadrupole/harmonic 
oscillator). Nevertheless, it has been shown that, for a more general Hamiltonian 
which includes our case as a special perturbative Dyson-like expansion for the 

Wigner function can be constructed which converges to the solution of the corresponding 
quantum Liouville equation. In general this approach, providing W{x,p,z) directly from 
the phase-space dynamics, should yield a smaller error. 

For a thin lens the Wigner transform, (|13]), of the approximate BWF, as given by (^) 
and (H), coincides with the approximate first-order solution of quantum Liouville equation 
0: the two approaches are therefore equivalent. As TWM has been mainly developed in 
the configuration space |^, and the approximate BWF for the typical potential used in 
particle accelerators has already been calculated (see, for example, |l|]-0), in this paper 
we have chosen to proceed with this latter approach. 

For what concerns the calculation of the higher-order solutions, solving directly ( [20| ) 
for W{x,p, z) should, in principle, be preferable to the procedure in which we first solve for 
BWF and then use (|13|), because less approximations are required, and therefore smaller 
errors are involved. Unfortunately it is very difficult to treat the Von Neumann equation, 
(PPP, numerically, especially if the classical potential contains powers in x higher than 
the quadratic one. In fact, in this case operators U x ± i{e/2)^ make (pO|) a partial 
differential equation of order higher than the second in the p-derivative, which is rather 
difficult to handle numerically. In order to treat BWF and Wigner transform numeri- 
cally, instead, one can profit of the very powerful methods developed for the Schrodinger 



equation, as it has been done, for instance, in Ref. [|T^. 

Although W is the distribution function of the system in the framework of TWM, 
due to well-known quantum mechanical properties, it is not positive definite. However, 
in QM it results to be positive for some special harmonic oscillator wave functions, called 



coherent states |Tl|, [T^, |T3|, , which give purely Gaussian density profiles. 

Coherent states for charged particle beams have been recently introduced in TWM 
in order to describe the coherent structures of charged particle distributions produced 



in an accelerating machine |T^. The fact that W can assume negative values for some 
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particular cases, reflects the quantum-like properties of both the wave function and the 
density operator. 

Remarkably, in the case of a pure quadrupole-like lens [U = kix'^/2) an interesting 
quantity, which estimates the r.m.s. area of the ensemble, can be expressed in terms of 
these momenta, namely 



71 



.1/2 



Since in this case the ground-like state is 

1 



IT 



X 



1/2 



exp 



the formula ( |T3D produces the following bi-Gaussian Wigner distribution 



(21) 



(22) 



Wq{x,p) = — exp 



2(T, 



(23) 







In addition, for the following non-coherent Gaussian-like state associated with the charged- 

exp +^ o.7./.A +^<P{z) , (24) 



particle beam |T5| 

^(x, z) 

where 



[2TTa\z)] 



1/4 



X 



X 



Aa^{z) 2eR{z 



kia 



dz^ ' 4a3 
definition (|13]) of W easily gives 







1 _lda 
R a dz 



dz 



4a2 



W{x,p,z) = — exp 



X 



2a^(z) 



X 



2cy\z) 



R{z) 



V 



(25) 



(26) 



Note that in this simple case of pure quadrupole the exact analytical solution (|23D and 
(p^) can be easily obtained directly by integrating (pO|). 

Note also that the argument of the exponential in (^61) is a quadratic form in the 
variables x and which can be written as 



F{x,p,z) = — ■y{z)x'^ + 2a{z)xp + l3{z)p'^ 
e L 



where 



a{z) 



aHz) 



zR{z) 



(3{z) 



oHz) 



7(z) 



+ 



aHz) 



Aa^z) eR\z 



These quantities are usually called Twiss parameters |TB|. By substituting 
(P^), we obtain a{z) = —{l/2)d(3{z)/dz, 7(2;) = (7p{z)/e and 



a^{z) = a{z) 



+ 



aHz) 



Aa'^iz) R^{z) Aa^{z) 



:i + 4a\z)) 



(27) 

(28) 
into (0)- 

• (29) 
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Consequently, we immediately get the following quantum-like version of the Lapostolle 
definition of emittance [|16] 



1 

{x^){p^) — -{xp + px)'^ = — = const, 
from which uncertainty relation of TWM can be easily derived 



Furthermore, the equation 



'y{z)x'^ + 2a{z)xp + l3{z)p^ 



(30) 



(31) 



(32) 



represents, for each z, an ellipse in the phase-space of area Tie/ 2 associated with the 
particle beam motion. 

Consequently, the operator 



J{x,p,z) = -f{z)x + 2a{z) + -f{z)p 



(33) 



is the quantum-like version of one of the well-known Courant-Snyder invariant |T^. In 
fact, it easy to prove that 

dj 



^ dz ~'~ 



J,H 







(34) 



where H is the Hamiltonian operator for the case of a quadrupole. 



3.2 Husimi function 

The phase-space description of a quantum system can be done in terms of another function 
which has been introduced by Husimi 0, the so-called Q-function or Husimi function 
(HF). In order to give an analogous definition of HF for charged particle beams in the 
context of TWM, we still replace Planck's constant with e, obtaining 

Q{xo,po,z) = - — / / Q*{u,z;xo,Po)p{u,v,z)Q{v,z;xo,po) du dv , (35) 

where Q{x, z; xo,po) is a coherent state associated with the charged particle beam defined 



as |T5 | 



e(u,z;xo,Po) 



[2nai] 



211/4 



exp 



{u - Xo{z)f 



with 



xoiz) = {u) 



Po{z) = {p) 



u |6|^ du , 



+ -po{z)u - i6o{z) 



(36) 

(37) 
(38) 
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and 



dSo 
dz 



2e V 14^2 



+ 



(39) 



Note that here Xq and po play the role of classical phase-space variables. 
By substituting (||) and ( ^61) in (|35|) we obtain 



Q{x,p, z) 



1 



27re-^27rcro 
X exp 



exp 



X 



2al 



dv 



du 



xiu + f ) i , 



(40) 



where for simplicity we have replaced the classical variables Xq and po with x and p, 
respectively. 

However, the integration over p and over x, does not give configuration-space and 
momentum-space distributions, respectively. But this function removes the pathology of 
negativity exhibited by W, and thus gives a more accurate description in the phase-space 
region where W is negative. 
Some properties of Q are in order. 



i) It is easy to prove the following normalization relation 

/OO /-OO 
/ Q{x,p, z) dx dp = 1 
-00 J — QO 



(41) 



ii) Definition (^Oj) can be cast in the following more convenient form 



1 



r 



Q{x,p,z) = ^= / dy exp 

Note that (H^) clearly shows that Q is positive definite 



I — W X H — ,z 
2t V 2' , 



(42) 



iii) Definition (|35|) of Husimi function, or equivalently (^), does not give in general a 
phase-space distribution coinciding with WF. But, in the case of a Gaussian BWF they 
must coincide, because in this case both Q and W reduce to the classical distribution. In 
particular, we observe that since in (^) or in (^) the constant ctq is involved, the present 
definitions of Q are suitable only to describe the phase-space distribution of eigenstates: 
this way Q does not explicitly depend on z. The natural generalization of Q for z- 
dependent BWF will be introduced later. Now we point out that, in connection with 
the Gaussian BWF given by (|^), Eq. (|42D gives a bi-Gaussian Q-function that does not 
coincide with (^) because of a scaling disagreement. This problem can be easily removed, 
by introducing the following new definition for Q 



Q{x,p,z) 



^7reA^/27^ag 



dy exp 



py 



2eA 



^ f VlT^x + ^, 



(43) 
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for any real number A. This is equivalent to introduce the following substitutions in the 
definition ( ^5]) 



X 



X 



p — — p 



(44) 



Under these substitutions, the normalization condition (^) is preserved. In order to 
symmetrize (^41), we choose A = 1. Consequently, 



X 



V2 



X 



P 



V2 



P 



(45) 



and 



becomes 



Q{x,p,z) 



ATie\l2'nol 



dy exp 



py 



^ i V2x + ^,z 



(46) 



Thus, by using (|^) in ( ^q) we obtain 

Q{x,p) = — exp 
7re 



X 



(47) 



which now coincides with the corresponding Wigner function (|23|). 

iv) For the non-coherent Gaussian-like BWF (|2^) , Eq. ( ^6|) is no longer valid, because it 
does not coincide with (^). In particular, at any z, the phase-space ellipses associated 
with W do not coincide with the contours of (|46|) for BWF given by (|2^) . This problem 
can be overcome by generalizing definition ( ^61) with the following 



Q{x,p,z) 



1 



dy exp 



16(T^{z) 8eR{z) V2e 



By substituting (|23) in (E3) we easily obtain 



py] ^{V2x+^,z 



(48) 



Q{x,p,z) = — exp 
Tie 



X 



2aHz) 



X 



2o\z) 



R{z) 



P 



(49) 



which now coincides with the corresponding WF given by 



(v) Note that the generalized definition (^81) suggests to start from (|35D where the 
coherent states 0(m, z; a;o,Po) are replaced by the following non-coherent Gaussian-like 
states 



eo(M,^;xo,Po) 



1 



\2t:g\z)\ 



1/4 



exp 



{u-x^{z)f % -R ( \ , (w- 2^0(2;))^ 



(50) 

This way if we start from ( |5BD with (|5DD, it is very easy to prove that it is equivalent to 
(^) for any Consequently, the generalized definition of ( pS]) with (|50D produce the 
same results of Wigner function definition (0) for any Gaussian-like states. 
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vi) Note that ( pOD are the non-coherent fundamental modes of the following set of solutions 
of (I) in the case oi U = {l/2)Kix'^ 



Qn{x,z;xo,Po) 



V2' 



Hn 



X — Xq[Z 



V2, 




a[z 



exp {i2n(j){z)) 



(51) 



These modes are analogous to squeezed and correlated Fock's states used in quantum 
optics [|l8l. In the next section we check if the equivalence between 



with (|50[) and 

(p^) is still valid for non-Gaussian BWF. In particular, we make a comparison between 
definition (0) and (|35D with (|50|) for the case of BWF given by @ and (|), and in addition 
we compare them with the corresponding results of a multi-particle tracking simulation. 



4 Numerical analysis 

A numerical study has been pursued in order to compare the description of the phase- 
space as given by TWM, both by means of WF and by means of HF, with the one resulting 
from standard particle tracking. 

A Gaussian flat (1-D) particle beam has been used as starting beam, with emittance 
e = 120 X lO-^m rad, and ctq = 0.05 m. From Eq. (^) it follows that, if a(0) = at the 
start, (Tpo = crp(O) = 1.2 x 10^'^ rad. (Note that the definition of emittance given in ([30| ) 
differs from the definition used in classical accelerator physics by a factor 1/2.) 

A simple device made of a quadrupole magnet plus a drift space has been considered 
as beam transport line: in addition sextupole and/or octupole aberrations have been 
included in the quadrupole. 

The Wigner function ([13D and the Q-function (|4^) have been computed by numerical 
integration for different combinations of aberration strengths, and have been compared 
with the results of the tracking of 7 x 10^ particles. 

Isodensity contours at 1, 2 and 3 a have been used to describe the particle distribution 
in phase-space, both before and after the passage through the simple device specified 
above. It is worth noting that with this choice only 2% of the particles are found beyond 
the contour at 2 a, and only 0.01% of them are beyond the contour at 3 a. 

In Figure 1. the starting distribution is shown together with the distribution emerging 
from the linear transport line. As already shown, in the linear case both WF and HF 
can be computed analytically (Eqs. ( ^61) and (^91) ) yielding exactly the same result of 
the conventional accelerator physics. Therefore, the output from tracking and the TWM 
predictions are in full agreement, and thus their superposition is not shown here. 

In the cases shown in Figure 2. non-linear perturbations have been added to the 
quadrupole lens: a sextupole perturbation in the first column, an octupole perturbation 
in the second column, and perturbations of both kinds in the third column. The results 
from the numerical computation of WF and HF are shown also superposed to the tracking 
results for a first set of perturbation values, corresponding to the distortions^ D = 0.125 
in the case of K2 = 0.06 and D = 0.042 in the case of K3 = 1.2 m~^. It should be 
noted that these values of distortion, chosen with the aim of enhancing the effects desired, 

^ The phase-space distortion can be defined as the ratio between the deflection Ap{x) = p{x) —p{xo), 
ai X = ao, and CpO: namely D = Ap{ao)/apo. It is easy to prove that D — 6t — 0^X2/ e for pure 
sextupole, and D ^ %lo ~ cTgA'a/Se for pure octupole. 
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are huge compared with what could be considered acceptable in any realistic single pass 
device. 

A quite good agreement with tracking can be observed for the contours at 1 and 2 o 
for both WF and HF. The contours at 3 o" show some discrepancies, more pronounced in 
the case of the Wigner function, which, as discussed in Section in the periphery of the 
distribution produces regions with negative phase-space density which yield an unrealistic 
distortion of the phase-space. By the use of the Q-function, instead, this effect is largely 
smoothed out, as shown in the last two raws of Figure 2.. 

In the plots displayed in the different columns of Figure 3. the same kinds of pertur- 
bations are included as in Figure 2., but with twice the strength, therefore yielding twice 
the distortion. 

These large values of distortion approach the limits of applicability of perturbation 
theory which, for the starting parameters and the quadrupole strength selected, are given 
by D < "ia^KxIe = 2.25 for the sextupole perturbation and D < Aao'^Ki/e = 3.00 for 
the octupole perturbation. 

Indeed larger discrepancies between TWM and tracking can be observed now, in par- 
ticular where the sextupole perturbation is present [D = 0.25). Nevertheless in these 
cases (columns 1 and 3 of Figure 3.) HF still describes fairly well the phase-space distri- 
bution for particles with amplitudes up to 1 — 1.5 a, whilst the WF distortion due to its 
negativity starts to show up already at amplitudes of the order of la. 

In spite of the discrepancies observed at large amplitudes due to the particularly strong 
perturbations used in this study, these results can be considered very satisfactory: the 
TWM phase-space description of beam dynamics in the presence of a non-linear lens, and 
in particular the one given by the HF, will be in more than reasonable agreement with 
the one given by classical accelerator physics for all realistic values of perturbation. 

5 Remarks and conclusions 

In this paper we have studied the phase-space distribution associated with the transverse 
motion of a charged particle beam travelling through a quadrupole-like device with sex- 
tupole and octupole deviations. In particular, a quantum-like phase-space analysis within 
the Thermal Wave Model has been developed and its predictions have been compared 
with the results of particle tracking simulations. 

To this end, we have first introduced the density matrix for the beam wave function 
(^. Then, following the usual definitions, we have constructed the Wigner transform 
(Eq. (|13D), and the Husimi transform (Eq. (p5D) with 9 given by (pOl). These functions 
represent the best candidates for describing the full phase-space evolution of the beam. 
Our aim was to compare the results of tracking simulations with the predictions of TWM, 
given in terms of W and Q, in order to enquire if the equivalence between (|35|) with (|50|) 
and (|T3|) is valid also for non- Gaussian BWF, and thus to determine the appropriate 
function to describe the phase-space dynamics. 

The Wigner transform of the BWF, which represents the natural quantum analogous 
of the classical phase-space distribution associated with the beam, has been numerically 
computed for the present problem. Its x and p-projections reproduce the configuration 
and momentum space distribution well , up to the first order of the perturbation theory. 
According to the results presented in section 4, this function reveals to be in good agree- 



12 



ment with the tracking results for small values of the integrated multipole-like strengths, 
for which conditions 0^X2! (3Ki) <^ 1 and alKzliyiKx) -C 1 are fully satisfied. For larger 
values of these parameters a little discrepancy appears. The reason is that in the last case 
the first-order perturbative expansion is not reliable enough. Nevertheless, we stress the 
fact that, beyond the contours at 2 o" and 3 a, only the 2% and 0.01% of the beam parti- 
cles are present, respectively. Summarizing, the comparison between the tracking results 
and the TWM predictions shows that: 

i) for small aberrations (conditions a^Kiji^K^ <C 1 and alK-ijiyiK]) <S 1 well satis- 
fied), both Q and W are in good agreement with the tracking results; 

ii) for larger aberrations, WF and HF exhibit the same order of discrepancy with respect 
to the tracking contours, but the distortions in the Wigner contours are much more evi- 
dent due to the negativity of this function which is responsible for a change in the contour 
concavity. 

In conclusion, up to the first order in the perturbation theory, and for thin lens ap- 
proximation, it results that Wigner function can be adopted as appropriate phase-space 
distribution in TWM. In fact it gives the correct x- and p-projections, which are in good 
agreement with the tracking configuration and momentum distribution, respectively. Fi- 
nally, although Husimi function does not give the right x and p-projections, according to 
our analysis it provides a better description as far as the only phase-space dynamics is 
concerned. 
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Figure Captions 



Figure 1. Phase-space distribution at the starting point (left) and after a quadrupole lens 
plus drift space (right). Starting from the center of the distributions the isodensity 
contours correspond to 1, 2 and 3 a, respectively. 

Figure 2. Comparison of TWM phase-space description with tracking (dotted lines) results. 

Parameter set No. 1. Starting from the center of the distributions the isodensity 
contours correspond to 1, 2 and 3 a, respectively. 

Figure 3. Comparison of TWM phase-space description with tracking (dotted lines) results. 

Parameter set No. 2. Starting from the center of the distributions the isodensity 
contours correspond to 1, 2 and 3 a, respectively. 
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